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ABSTRACT 

We use the Millennium Simulation to measure the cross-correlation between 
halo centres and mass (or equivalently the average density profiles of dark 
haloes) in a ACDM cosmology. We present results for radii in the range 
10 /i _1 kpc < r < 30 ft, _1 Mpc for halo masses in the range 4 x 10 10 h~ 1 M Q < 
M200 < 4 x 10 14 /i _1 M Q . Both at z = and at z — 0.76 these cross-correlations 
are surprisingly well fit by approximating the inner region by a density profile 
of NFW or Einasto form, the outer region by a biased version of the lin- 
ear mass autocorrelation function, and by adopting the maximum of the two 
where they are comparable. We use a simulation of the formation of galaxies 
within the Millennium Simulation to explore how these results are reflected 
in cross-correlations between galaxies and mass. These are directly observ- 
able through galaxy-galaxy lensing. Here also we find that simple models can 
represent the simulation results remarkably well, typically to <J 10%. Such 
models can be used to extend our results to other redshifts, to cosmologies 
with other parameters, and to other assumptions about how galaxies populate 
dark haloes. The characteristic features predicted in the galaxy-galaxy lensing 
signal should provide a strong test of the ACDM cosmology as well as a route 
to understanding how galaxies form within it. 
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1 INTRODUCTION 

Weak gravitational lensing has opened a new window 
onto the large scale distribution of matter. Gravita- 
tional lensing by foreground mass induces correlated 
distortions, or shear, in the observed shapes of distant 
galaxies. In galaxy-galaxy lensing, the signal from many 
galaxies is added together in order to measure the av- 
erage (projected) distribution of mass around galaxies. 
This can be interpreted as the mass in the extended 
dark matter halos which surround galaxies, or, more 
generally, as the cross-correlation between lens galax- 
ies and the projected mass distribution. Several groups 
have successfully applied this technique to large imag- 
ing surveys to derive constraints on the mass associated 
with galaxies as a function of g alaxy propertie s such 
as luminosity and mor p holog y dBrainerd et al.l 
dell' Antonio and Tvsonl Il996l : iGriffiths et al ' 



1996 



1996 



Hudson et al.l Il998l; iFischer et al.l l200d: iMcKav et al l 
20011 ; iHoekstra et alj|2003l ; ISheldon et al.|[20ol) . 

Theoretical predictions for cross-correlations be- 
tween galaxies and mass, £ gm , have made use both of nu- 



merical simulations (iGuzik and Selia kll200ll;lYang et al.l 
120031 ; iTasitsiomi et alj|2004 IWeinberg et alj|2004) and 
of an a lytic halo models (ISeliakll2000l ; IGuzik and Selialj 
120021 ). ITasitsiomi et all (|2004n show that the amplitude 
and shape of £ gm predicted by cosmological simulations, 
when combined with a simple model for populating 
halos with galaxies, are in good agreement with ob- 
servational re sults based on Sloan Digital Sky Survey 
(SPSS ) data (|Sheldon et all 120041 1. iMandelbaum et~al1 
|2005l ) also find that these simulation results can be ac- 
cur ately reproduced by the h alo model of ISeljak (2000) 
and lGuzik and Seliakl (|2002f) . 

In this work we calculate the cross-correlation be- 
tween halos and mass, £hm, and between galaxies and 
mass, £ gm , in the Millennium Simulation, a very large, 
high-resolution simulation of a ACDM cosmology. We 
also present simple models for £hm and £ gm which can 
be used to interpret the shapes of these functions and 
to extend our results to other redshifts, cosmologies and 
halo population models. 

The organization of this paper is as follows. In Sj2] 
we describe the Millennium Simulation and the halo and 
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galaxy catalogues used in this study. In Sj3]we show £hm 
calculated for halo samples with a range of masses and 
we present a model which accurately reproduces these 
results. In Sj4]we show £ gm for both central and satellite 
galaxies as a function of their luminosity, and we present 
models for each of these and for the combined galaxy 
sample. Finally, we summarize our results in [J5] 



2 SIMULATIONS 

This study makes u se of the Millennium Simulation^ 
(|Springel et al.ll2005l ). a large cosmological N-body sim- 
ulation carried out by the Virgo Consortiurr0. In this 
simulation a flat ACDM cosmology is adopted, with 
f2dm = 0.205, Qh = 0.045 for the current densities in 
cold dark matter and baryons, h = 0.73 for the present 
dimensionless value of the Hubble constant, as = 0.9 
for the rms linear mass fluctuation in a sphere of radius 
8 /i _1 Mpc extrapolated to z = 0, and n = 1 for the slope 
of the primordial fluctuation spectrum. The simulation 
follows 2160 3 dark matter particles from z = 127 to z = 
within a periodic cube I/box = 500 /i _1 Mpc on a side. 
The individual particle mass is thus 8.6 x 10 s /i _1 M Q , 
and the gravitational force is softened with a Plummer- 
equivalent comoving softening of 5 h~ 1 kpc. Initial con- 
ditions were generated using the Bol tzmann code CMB- 
FAST ijSeliak and Zaldarriagalll996l ) to generate a real- 
ization of the desired power spectrum which w as then 
impo sed on a glass-like uniform particle load (jWhitei 
1 19961 ). A modifie d version of the TREE-PM N-body 
code GADGET2 dSpringel et al.ll2001bl. [iooEh was used 
to carry out the simulation and full particle data are 
stored at 64 output times approximately equally spaced 
in the logarithm of the expansion factor at early times 
and at roughly 300 Myr intervals after z = 2. 

In each output of the simulation, halos are iden- 
tified using a friends-of-friend s (FoF) groupfind er with 
a linking length of 6 = 0.2 ijPavis et alJ Il985h . Each 
FoF halo is decomposed into a collection of locally over- 
dense, self-bound subst ructures (or subhalos) using the 
SUBFIND algorithm of lSpringel et alJ (|2001ah . Of these 
subhalos, one is typically much larger than the others 
and contains most of the mass of the halo. We identify 
this as the main subhalo and define its centre as the po- 
sition of the particle with the minimum potential. The 
virial radius, V2oo, is defined as the radius of a sphere 
that encompasses a mean density 200 times the critical 
value, and the virial mass, A/200, is the mass within this 
radius. 

Semi-analytic techniques have been used to simu- 
late the evolution of the galaxy populat ion within the 
Millen nium Simulation, as describ ed by ISpringel et al.l 
(120051 ) and ICroton et ah! (|2006l ). In this approach 
the evolution of the baryonic component is followed 
using a set of simple prescriptions for gas cooling, 
star formation, supernova and AGN feedback, chem- 
ical enrichment, galaxy merging and other relevant 



1 http:/ /www. mpa-garching.mpg.de/millcmiium/ 

2 http://www.virgo.dur.ac.uk/ 



physical processes. These models can be applied 
repeatedly to the stored histories of the dark matter 
halos and subhalos with different parameter choices 
for the model, or indeed with different physical as- 
sumptions in the modelling. In this work we use the 
semi-analytic galaxy catalogue described in detail in 
ICroton et al.l |2006). We note that various aspects 
of the clustering properties of halos and galaxies in 
the Millennium Simu l ation have been in vestigated in 
ISpringel et~ai1 (120051) iGao et alJ (I2005I). lHarker et alJ 



ICroton et al 



(2007) and 



(12009) I Wang et alJ J2006l). ISpringel et alJ (120061) 



(l2006l). iLi et all <|2006| ), ICroton et al l 
Gao and White! (|2007i ). Studies of halo 



de nsity profile s in th e sim ulation have been carr ied out 
bv lNeto et all (120071 ) and lGao and White! (|2007l ). 



3 HALO-MASS CROSS-CORRELATIONS 

Given a density field, p(x), the density fluctuation field 
is defined as 



<5(x) = 



P(x) - p 



(1) 



The two-point autocorrelation function is defined as 

f(r) = <<5(x)5(x + r)). (2) 

Assuming isotropy reduces this to a function of the sep- 
aration, £(r). This is interpreted as a measurement of 
the excess probability above random of fi nding a pair 
of objects with separation r (jPeeblesll 19801 ). In the case 
of two different populations of objects, <5 a (x) and <5b(x), 
the two-point cross-correlation function is given by 



fab(r) = (5 a (x)(5b(x + !•)). 



(3) 



If we consider the cross-correlation between halo centres 
and mass, the cross-correlation function £hm(r) simply 
reflects the spherically-averaged halo density profile av- 
eraged over all halos in the sample. This can be seen by 
combining eqs. [1] and [3] to give 

(p(r)> - p m 



(4) 



where r is the radial distance from the halo centre and 
Pm = pcritfi m is the mean density of the Universe. 

Figure [T] shows halo-mass cross-correlations along 
with the mass autocorrelation function for the Millen- 
nium Simulation. The halo-mass cross-correlations are 
computed using the centres of the main subhalos of all 
FoF halos with mass M200 > 4 x 10 10 ft _1 M , cor- 
responding to J> 50 particles. Seven halo samples are 
shown, each spanning a factor of two in mass. The num- 
ber of halos in each sample is listed in Table [1] We have 
computed cross-correlations also for the mass ranges not 
included in this plot and table, and we have checked that 
they also fit the models we discuss below. However, for 
clarity we refrain from showing them or discussing them 
further. 

The shape of the halo-mass cross-correlation is 
clearly defined by two parts, commonly referred to as 
the one-halo and two-halo terms since they are domi- 
nated by particles within the same halo and in different 
halos, respectively. Figure [1] shows that on large scales 
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Table 1. Properties of halo samples taken from the Millennium Simulation. 



Mass range 


-*» halos 


< A/200 > 


< #200 > 


Fit M 200 


Fit R200 


Fit C200 


''"trans 






[/i- 1 M ] 


h _1 Mpc 


[h- 1 M ] 


ft, _1 Mpc 




h _1 Mpc 


[4.0,8.0] x 10 10 


2491565 


5.59 x 10 1(1 


0.062 


5.66 x 10 10 


0.062 


9.8 


0.11 


[1.6,3.2] x 10 11 


774837 


2.23 x 10 11 


0.098 


2.26 x 10 11 


0.099 


9.1 


0.19 


[6.4, 12.8] x 10 11 ] 


230429 


8.92 x 10 11 


0.16 


9.05 x 10 11 


0.16 


8.3 


0.33 


[2.56,5.12] x 10 12 


67400 


3.56 x 10 12 


0.25 


3.62 x 10 12 


0.25 


7.5 


0.74 


[1.02,2.048] x 10 13 


18322 


1.42 x 10 13 


0.39 


1.45 x 10 13 


0.40 


6.7 


1.28 


[4.1,8.192] x 10 13 


4274 


5.63 x 10 13 


0.62 


5.79 x 10 13 


0.63 


5.7 


2.20 


[1.64,3.277] x 10 14 


650 


2.22 x 10 14 


0.98 


2.32 x 10 14 


1.00 


4.7 


3.78 



10° 
10 5 
10" 
10 3 
10 2 
10' 
10° 

10-' 
10- 2 




, 4.0x10'°<M 200 <8.0x10 
1.6x10"<M 200 <3.2x10 v 
6.4x1 0"<M 200 <1.2x10 ,: 
2.56x10 ,2 <M 200 <5.12x10' 2 
1. 024x1 ,3 <M 200 <2.048x10' 3 

, 4.096x10 ,3 <M 200 <8.192x10' 3 

1 ,6384x 1 o"<u m) <3.2768x 1 0" 
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r (Mpc/h) 
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Figure 1. Cross-correlations between halo centres and the 
mass in the Millennium Simulation at z = 0. The solid curve 
shows the mass autocorrelation function, § m m for compari- 
son. On large scales, r > 3 h _1 Mpc, the cross-correlation 
functions follow the mass autocorrelation function with a 
constant bias factor. On small scales, fh m (r) reflects the 
mean density profile of the individual halos. Note the sharp 
transition between the two regimes. 



Chm(^) follows closely the mass auto-correlation func- 
tion, £ mm , with a mass-dependent offset in amplitude 
known as the halo bias factor, b(M). The transition be- 
tween the two regimes is remarkably sharp and takes 
place at an overdensity of approximately seven times 
the mean density, i.e., £hm — 6. 

We construct a simple model for (i, m (r) as follows: 



Cmodel (r;M) 



Cih(r) if£ihM>6h(r), 
6t(r) if £ih(r) < 6h(r), 
Phaio(r;M) - p m 



(5) 

(6) 
(7) 



where £ijn(r) is the mass autocorrelation function pre- 
dicted by linear theory. The main ingredients of the 
model are the halo density profile, phaio, and the bias 
factor 6(M) which we now describe in detail. 



3.1 The one- halo term 

The density profiles of CDM halos have been stud- 
ied extensively with high resolution N-body simula- 
tions over the past decade. Early results indicated 
that the densi ty increases steep ly towards the cen- 
tres of halos dFrenk et all Il985l; iQuinn et all | l98£ ; 

23D- " 



Dubinski and Carlberd Il991f ). iNavarro et all 1(1996 
19971 . hereafter NFW) suggested the following simple 
fitting formula to describe the density profile of simu- 
lated halos: 

PNFw(r) _ So , . 

Pcrit (r/r s )(l + r/r s ) 2 ' { > 

where p cr it is the critical density. Note that the slope of 
the NFW profile is shallower (steeper) than the isother- 
mal profile inside (outside) the characteristic scale ra- 
dius, r s . Integrating this density profile out to the virial 
radius, r2oo, gives the following relation for the dimen- 
sionless density parameter 



8 



200 



c 200 



(9) 



3 ln(l + C200) — C2oo/(l + C200. 

where the concentration parameter C200 = T20o/r 3 . Since 
the halo mass and virial radius are related through 
-M200 = 200 pcrit (4-7r/3) r| 00 , the independent param- 
eters in the NFW profile are effectively the halo mass 
and concentration. Furthermore, these properties are 
known to be correlated in simulated halos, in the sense 
that low-mass halos are more central concentrated than 
high-mass halos. This is generally interpreted in terms 
of the mean density of the universe at the time of for- 
mation of a halo. Since low-mass systems typically col- 
lapse at higher redshift, the characteristic density and 
concentration of such systems is larger with respect to 
high-mass systems. The concentration-mass relation has 
been studied extensively with cosmological simulations 
and several authors have proposed models for predict- 
ing the avera ge value of C200 as a function of halo mass 
and redshift (INavarro et al.lll997l: iBullock et aD | 200ll : 
Eke et all l200ll: iMaccio et all 120071; iNeto et"afl 12OO7I; 



Gao and Whitej[2007l ) . Adopting such a model therefore 
fully specifies the density profile of a typical halo of a 
given mass M200 at a given redshift. 

The most recent high resolution N-body simula- 
tions have revealed small but significant deviations from 
the NFW formula. INavarro et ail (|2004 hereafter N04) 
show that the density profiles become shallower towards 
the halo centre more gradually than the NFW formula 
predicts, causing NFW fits to underestimate the density 



© 0000 RAS, MNRAS 000, 000-000 



4 E. Hayashi and S.D.M. White 



in the inner regions. These authors propose an improved 
fitting formula with the following form: 



ln^ = ^ 

P-i ct 



r-2 



- 1 



(10) 



where p_2 and r_2 are the density and radius at 
which the logarithmic slope of the density profile 
dlogp a /dlogr = —2, and the parameter a, controls 
the rate of change of the logarithmic slope with ra- 
dius. Higher values of a cause the profile to become 
shallower more quickly toward the centre of the halo. 
Unlike the NFW profile, eq. \W\ does not converge to 
a power law at small radii, instead reaching a finite 
density at the centre. N04 found a relatively small 
range of values, 0.12 < a < 0.22, to fit the density 
profiles of simulated halos over a wide range in mass 
(8 x 10 9 /i _1 M <; M 200 <; 8 x 10 14 /i _1 M ). A spatial 
density profile of this fo r m wa s first proposed in a differ- 
ent context bv lEinastd \l96$j so we will hereafter refer 
to it as the Einasto profile. Prada et al.l (|2006l ) recently 
used the NFW and Einasto formulae to fit halo density 
profiles out to radii far beyond the virial radius. These 
authors find the mean density to give a significant con- 
tribution to the density profile at large distances, i.e., 
the halo profile is better fit by 

phalo(r) = p a (r) + p m . (11) 

With this modification, the one-halo term becomes 

P«(r) 



fih(r) = 



(12) 



In subsections 13.31 and 13.51 we investigate the accuracy 
of our model for £hm using the NFW and Einasto fitting 
formulae in the one-halo term of the model. 



3.2 The two-halo term 

The two-halo term of the model is specified by the 
bias factor, b(M), and the mass autocorrelation func- 
tion calculated from linear theory, £ii n (r). Halo bias has 
been studied extensively in the context of hierarchi- 
cal structure formation scenarios. Assuming a Gaussian 
distribution of initial density fluctuations specified by 
a(M), the rms linear mass fluctuation (extrapolated to 
z = 0) within spheres that on average contain mass M, 
iMo and~W hitc ( 1996) derive an analytic model for halo 
bias, b(v, z), where 



S c (z) 



a(M) 



(13) 



is the dimensionless amplitude of fluctuations, or peak 
height, that produces halos of mass M at redshift z and 
5 c (z) is the linear overdensity (again extrapolated to 
z = 0) for which a spherical perturbation would collapse 
at redshift z.. The characteristic halo mass for clustering 
M*(z) is defined by cr(M«) = 8 c (z) and halos more (less) 
massive than M* are more (less) strongly clustered than 
the underlying mass density field. For the Millennium 
Sim ulation cosmolo g y, M t ( z = 0) = 6.15xl0 12 h~ l M & . 

IMo and White! l|l996h showed that their model 
accurately describes the bias in the autocorrelation 
function of dark matter halos with respect to the 



that of the mass in cosmological N-body simulations. 
Further testing against higher resolution simulations 



(Jine 1998; iGovernato et al.l Il999l; ISheth and Tormen 


1999 


; Kravtsov and KlvDin 19991; Colbere et al. 


2000 


: ISheth et al 


l200ll: ISeliak and Warrenl 20041: 


Mandclbaum et al. 


l2005h. Most recently, iGao et al. 


(200 


5|) compared 


these models with halo bias in 



the Millennium Run. Figure 1 of IGao et al] i|2005l) 
shows that the Millennium Run res ults are reasonably 
well-m atched by the bias model of ISheth and Tormenl 
(jl999l ) 



b{v) = 1 + 



+ 



2p 



(14) 



Soil + {av 2 )P) ' 

with the paramete r value s a = 0.73 and p = 0.15 of 
iMandelbaum et alj (l2005h b ased on fits to the simula- 
tions of ISeliak and Warrenl (|2004T ). We therefore adopt 
eq. [T?] with these parameter values as the bias formula 
in the two-halo term of our model for £hm- 



3.3 Model fitting results 

Figure [2] shows the accuracy of our model when we 
choose the halo density profile to be of NFW form. The 
downward pointing arrows indicate the transition scale 
^trans between the one-halo and two-halo regimes of the 
model, also listed in Table [T] We take the halo mass 
in the model to be the geometric mean of the upper 
and lower halo mass limits in each halo sample, i.e., 
A^modei = (MioMhi) 1 / 2 . Here we take the concentration 
to be a free parameter and compare our best fit values 
to the concentration-mass relations proposed by various 
authors. The best fit value is obtained by minimizing the 
root mean square of (£ hm — £modd)/£hm. 

Several trends are apparent in the accuracy of our 
model fits. In the one-halo regime, the NFW profile fits 
the halo density profile to within about 10%. There is 
some indication of the systematic "U-shaped" residuals 
found by N04 for fits to individual halo density profiles 
with the NFW formula. This indicates that the shape 
of the NFW profile does not perfectly capture the sim- 
ulation results. We return to this issue later, when we 
examine the accuracy of fits using eq. [10] to model the 
halo density profile. 

In the two-halo regime, the deviations in the model 
at large separations, r J> 3 /i~ 1 Mpc, are dominated by 
the quasi-linear distortion due to the large-scale move- 
ment of halos with respect to each other. This is illus- 
trated by the solid line in Figure [2] which shows the ratio 
of the mass autocorrelation function to the linear the- 
ory prediction, £mm/£iin- This distortion has been inves- 
tigate d in numerous studies of the matter power spec- 
trum fecliak 200^ lMaand Frvll200d ; [Scoccimarro et all 



120011 ; ISmith et all 120031 ; ICole et all 120051 ) and various 
correction factors have been proposed. In this work we 
prefer to neglect the distortion an d adopt the simpler 
assum ption of pure linear theory (|Peacock and Smith! 
2000). This approach simplifies the calculation of the 
model and avoids the introduction of the additional pa- 
ra meters that are pre sent i n correction facto rs proposed 
bv lSmith etafl (|2003T ) and lCole et"afl (120051 ) at the cost 



© 0000 RAS, MNRAS 000, 000-000 



Halo-mass and galaxy-mass cross- correlations 5 



E 
E 



0.3 



0.2 



0.1 



-0.0 



-0.1 



-0.2 



-0.3 



1 



f hrr 
£hir 

?hrr 

£hnri 




4.0x10'?<Mj^ <8.0x10" > 

i.6xio"<M 20 o<3.2xio" 

6.4x1 0"<M 200 <1.2x10' J 
2.56x10 l2 <M 200 <5.12x10' 2 
1. 024x1 ,3 <M JOO <2.048x10' 3 
4.096x1 ,3 <M 200 <8.1 92x1 0' 3 
1. 6384x1 0"'<M 200 <3.2768x10" 
' ' 



0.01 



0.10 1.00 
r (Mpc/h) 



10.00 



Figure 2. Deviations between our measured halo-mass cross- 
correlations and the simple model given by eqs. I5I7I as- 
suming the NFW profile in the one-halo term. On large 
scales, r > 3 /i _1 Mpc, the deviations are dominated by a 
quasi-linear distortion that is apparent in the ratio of the 
mass autocorrelation function to the linear theory prediction, 
£mm/£ii n (solid line). For low mass halos, M200 < 3.2 X 10 11 , 
the model also fails on intermediate scales suggesting a cut- 
off in power should be applied to the linear prediction on 
these scales. Otherwise the accuracy of the model better than 
~ 10%. 



of a systematic residual in our model fits on the order 
of 25% at ~ 10 /i -1 Mpc. 

In the two lowest mass halo samples shown in Fig- 
ure [21 the two-halo term in the model appears to over- 
predict £hm at overdensities corresponding to £hm > 7. 
This suggests that the linear theory prediction should 
also be modifie d by a cutoff in p ower at small scales, 
as suggested bv lSmithet al] l|2003l ). We investigate this 
in further detail in Figure [3] where we exa mine the lin- 
ear an d quasi-linear terms predicted by the lSmith et al.l 
l|2003h HALOFIT model for the Millennium Simulation 
cosmology. The correlation function is related to the 
power spectrum by the integral relation 

. , , f . 9 , . . sin kv elk 
£(r)= / A 2 (k)—— — , (15) 
J kr k 

and the HALOFIT model power spectrum is given by 
the sum of the quasi-linear and one-halo terms, 



A 2 (fc) = A|(fc) + A|(fc), 



(16) 



where the quasi- linear term Aq(fc) includes an ex- 
ponential cutoff at the nonlinear wavenumber fc CT ~ 
0.3 h Mpc -1 for ACDM . 

The top panel of Figure [3] shows that the correlation 
function calculated from the HALOFIT model power 
spectrum provides a good fit to the mass autocorrela- 
tion function £ m m- The cutoff in power at high k built 
into the HALOFIT model corresponds to a flattening 
of the corresponding correlation function, £q with re- 
spect to the linear theory prediction, £n n at an overden- 



sity £q ~ 4. ISmith et al.l (|2003h provide fitting formu- 
lae for Aq(fc) and Ajj(fe), each with numerous param- 
eters tuned to give a good match to the matter power 
spectrum in various CDM-based cosmological N-body 
simulations. However, Figure [3] shows that £q cannot 
be used to improve the model for £hm of low-mass ha- 
los. This is illustrated in the bottom panel, where we 
show the ratios £hm/£Q and £hm/£im for halos with mass 
4.0 x 10 10 h^Mo < M200 < 8.0 x 10 10 h~ 1 M Q . We find 
that using the quasi-linear term from the HALOFIT 
model actually produces a marginally poorer fit to the 
measured shape of halo-mass cross-correlations. The 
HALOFIT model imposes a cutoff in power at the 
non -linear wave n umbe r k a . Our results suggest that 
the ISmith et al.l (120031 ) prescription adopts a value of 
k a that is too low in the present context. Modifying 
their prescription for the non-linear scale would require 
the recalibration of all of the fitting formulae in the 
HALOFIT model in order to recover the fit to the mass 
autocorrelation function. This is beyond the scope of 
the present paper. 



3.4 Halo concentrations 

Figure [4] shows the best-fit NFW concentration param- 
eter values versus halo mass for our model fits. Also 
shown are the pred ictions of t he co ncentration-mass 
models propo s ed by lEke et al.l l|200ll . hereafter ENS), 
iBullock et all (l200ll. here a fter B 01), and the power- 
law fits of Maccio et all l|2007l . hereafter M07) and 
iNeto et al. I (|2007l ). Our best-fit values appear to follow 
the B01 and M07 results for halo masses M vir < 5 x 10 12 , 
but at higher masses the results are better described by 
the ENS model. We find that a power-law provides a 
reasonable fit to the data, with the same normaliza- 



tion, but a slightly shallower sl o pe (c oc M 080 ) than 
found by M07 and INeto etafl (|2007r i (c oc M -0109 ). 
We note that the difference between all these models 
is relatively small with respect to the 1 — a scatter, 
Alog(c v i r ) = 0.14, found for individual halos in B01. 
Since the ENS model provides a reasonably close match 
to our best-fit concentration values and has been cali- 
brated for redshifts z > and for different cosmologies, 
we adopt it as our halo concentration model hereafter. 



3.5 Halo density profile 

In this subsection we investigate whether our model is 
improved if we change the choice of halo density profile. 
We here use the Einasto profile (eg. I10p . the improved 
fitting formula proposed by N04. This contains an ad- 
ditional shape parameter, a, which controls the rate at 
which the slope of the density profile changes with ra- 
dius, with typical value s lying in the range [0.12, 0.22] 
(N04. |Prada et al.ll2006l ). 

Figure[S]shows the deviations of the measured halo- 
mass cross-correlations from our best fits using the 
Einasto formula. Here we focus on halo samples with 
M200 > 6.4 x 10 11 h^Mo, since deviations from the 
fits are then not dominated by the inadequacy of the 
two-halo term on intermediate scales (see Section f3.3|l . 
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Fi gure 3. Top panel: Mass autocorrelation function, £mm, 
and halo- mass cross-correlation functions, £hm> f° r low-mass 
halos (4.0 X 10 10 < M 20 o < 8.0 X 10 10 ). Also shown are the 
predictions of linear theory £u n and of the lSmith et al.l l|2003f ) 
HALOFIT model. Neither f Un nor the HALOFIT quasilinear 
model £q provide a good match to the shape of £hm on large 
scales. Bottom panel: The ratio between £h m and the linear 
and quasilinear predictions shows that the shape of fhm is 
marginally better fit by £u n than by £q . 



Figure 4. NFW concentration versus halo mass from £h m 
model fits compared with the predictions of various models. 
The solid line shows the best-fit power law with c oc M~ o oao . 
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Figure 5 . Deviations between our measured halo-mass cross- 
correlations and the simple model given by eqs. 1 5171 with 
the Einasto profile in the one-halo term. On small scales the 
Einasto profile provides a better fit than the NFW profile, 
with deviations < 5%. 



The deviations are reduced significantly by using the 
Einasto model, <C 5% compared to <C 10% for fits 
using the NFW profile. The best-fit values of a are 
shown as a function of halo mass in Figure [6] We find 
that a tends to increase with halo mass, ranging from 
a ~ 0.12 for M 20 o ^ 10 12 /i _1 M to a ~ 0.2 for 
M 20 o =^3x 10 14 ft _1 Mo. 



3.6 Higher redshift results 

Since galaxy-galaxy lensing studies are sensitive to 
lenses at redshifts well above zero, it is interesting to 
check our models against cross-correlations at z > 0. 
In the COSMOS survey, for example, the galaxy-galaxy 
lensing signal is significant for lenses over the redshift 
range 0.2 < z < 1.2 with the main contribution from 
z ~ 0.5 — 0.9 (A. Leauthaud, private communication). 
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Figure 6. Best-fit value of a from the Einasto fitting formula 
versus halo mass. Higher values of a correspond to density 
profiles that becomes shallower more quickly towards halo 
centre. The best-fit value of a tends to increase with halo 



Figure [7] shows £hm f° r the z = 0.755 output of 
the Millennium Simulation. The results are qualitatively 
similar to those at z = 0, and the bottom panel of the 
figure shows the deviations from the best fits using our 
(Einasto) model. We find that the deviations are on the 
order ±10% for the higher redshift £hm- Figure [6] also 
shows the best fit values of a for this redshift. As in 
the z — case, we find that a tends to increase with 
mass for Af 2 oo J£ 3 x 10 11 h~ 1 M Q - In the left panel 
of this figure we plot the best fit a values against the 
peak height defined by eq. 1131 Plotted in this way, we 
find that the 2 = and z = 0.755 results are in good 
agreement with each other. We note that a similar result 
has been obtained by Gao et al (in preparation) in a 
study of cluster halos from the Millennium Simulation. 
This provides a simple way to estimate appropriate a 
values for other redshifts and masses. 



3.7 Shear 

Galaxy-galaxy lensing studies measure the average tan- 
gential distortion or shear of background source galaxies 
due to foreground lenses. If the redshift distributions of 
the sources and lenses are known, the average tangen- 
tial shear, yt, is related to the surface mass density as 
follows 



E crit7t = E(< R) - E(i?) = AE(i?), 



(17) 



where E(< R) is the mean surface density within the 
projected radius R, and the critical surface density is 
given by 



4ttG D,D, t 



(18) 



where D s and Di are the angular diameter distances to 
the lens and the source, respectively, and Di s is the an- 
gular diameter distance between the lens and the source. 
In order to compare our results with observational mea- 
surements, we calculate the surface mass density by pro- 
jecting the three-dimensional mass density: 

Ap(r) 



E(R) = 2 f 

J r 

where Ap(r) = p(r) — p. 



R 2 



=rdr, 



(19) 
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Figure 7. Top panel: Halo-mass cross-correlations calculated 
at z = 0.755. Bottom panel: Deviations between J, m and our 
model assuming the ENS halo concentration relation and the 
Einasto density profile. Deviations from the model are on the 
order of 10 - 20%. 



Figure[8]shows the quantity AE(ii) calculated from 
our £hm curves and our model fits. Note that the model 
fits have not been recalculated to minimize the de- 
viations in AE(i?). We find that down to a scale of 
~ 0.03 - 0.05 /i _1 Mpc, below which our calculation of 
E(-R) from the simulation becomes unreliable, the de- 
viations from the model fits are <; 10% for both z = 
and z — 0.755. 



4 GALAXY-MASS 

CROSS-CORRELATIONS 

We now investigate cross-correlations between galaxies 
and mass in the Millennium Simulation. As described 
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dynamical friction time, after which they merge with 
the central galaxy of their halo. They are a significant 
fraction of the fain t galax ies i n the mode l . iGao et al] 
(|2004MWang et all (120061 1 and lSales et~aH i|2007l ) show 
that they must be included if the model is to predict 
the small-scale clustering of galaxies accurately. 



4.1 Cross-correlations between central 
galaxies and mass 

We first consider cross-correlations between central 
galaxies and mass. Since central galaxies are placed on 
the most bound particle of their host halos, £ gm for cen- 
tral galaxies is equivalent to £hm averaged over the corre- 
sponding sample of FoF halos, and so should be similar 
to the £hm presented in Section [3] provided that central 
galaxies are selected according to a property that cor- 
relates well with host halo mass. Figure [9] shows £ gm 
for central galaxies selected according to r-band abso- 
lute magnitude in the range —20 < M r < —24. These 
cross-correlations indeed appear very similar in shape 
to those of Figure [T] 

We apply the simple model presented in Section [3] 
based on the Einasto profile of eq. 1101 We take the mass 
of the host halo as a free parameter which determines 
the bias factor according to eq. [14] and the radius and 
concentration of the mean halo density profile according 
to the ENS model. The best fit value of the halo density 
profile parameter a ~ 0.15 for the central galaxy £ gm . 
This is quite similar to the values found above for halos 
in the relevant mass ranges. 

The distribution of host halo masses for central 
galaxies in the luminosity range —21 < M r < —20 is 
shown in Figure 1101 This distribution is highly asym- 
metric and ranges over four orders of magnitude in mass. 
Nevertheless, the best fit model recovers the mean halo 
mass within about 30%. The bottom panel of the figure 
shows mean and best-fit halo mass values as a function 
of M r . The fit values tend to underestimate the mean 
halo mass, typically by about 30%. 



Figure 8. Top panels: AS(iJ) calculated from £i lm at z = 
and 2 = 0.755. Solid lines show AS(H) corresponding to 
best fits to §hm- Bottom panels: Deviations from our fits. The 
downturn at small scales r < 0.03 h~ 1 Mpc is due to integrat- 
ing over scales below the resolution limit of the simulation. 
On scales larger than this, the magnitude of the deviations 
are on the order of 10% in terms of AE(il). 



in Section [5] semi-analytic models have been grafted 
onto the simulation in order to predict the evolution of 
the galaxy population. A central galaxy is associated 
with the main subhalo of each FoF halo and satellite 
galaxies are associated with the other subhalos. In ad- 
dition, there are so-called 'orphan' galaxies. These are 
satellite galaxies whose associated subhalos have been 
tidally stripped to below the resolution limit of the sim- 
ulation. Such a galaxy remains identified with the in- 
dividual particle which was the most bound member of 
the subhalo at the last time it could be identified. Or- 
phan galaxies are assumed to survive for an estimated 



4.2 Cross-correlations between satellite 
galaxies and mass 

We now consider cross-correlations between satellite 
galaxies and the mass. These galaxies were once asso- 
ciated with the dominant subhalo of a FoF group but 
are now centred on a secondary subhalo within a larger 
FoF group and may have experienced significant mass 
loss due to tidal stripping. Indeed, as noted above, a 
significant fraction are 'orphans' and have lost their 
subhalo entirely, remaining associated with the parti- 
cle which lay at subhalo centre when it was last identi- 
fied. Observation shows that the baryonic components 
of galaxies are substantially overdense with respect to 
the dark matter, and simulations suggest th at they are 
therefore m ore resistant to tidal disruption (Kat z et al.l 
Il992l .[l996). The semi-analytic galaxy formation model 
takes this into account by allowing orphans to survive 
for a dynamical friction time before merging them with 
the central galaxies of their halos. 
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Figure 9. Top panel: Galaxy-mass cross-correlations at 
z = for central galaxy samples selected according to r- 
band absolute magnitude. Bottom panel: Deviations from our 
model for halo-mass cross-correlations, where the mean host 
halo mass has been taken as a free parameter. 



Table [2] shows the numbers of central and satel- 
lite galaxies i n the Millennium Simulation model of 
ICroton et all (2006) as a function of M r . Although 
there are fewer satellites than central galaxies for all 
Mr < —17, the fraction of satellites increases with de- 
creasing luminosity, and is about 40% of the galaxy 
population for —19 < M r < —17. The fraction of or- 
phans also increases with decreasing luminosity, and is 
more than half of the total satellite population in the 
— 18 < Al r < — 17 magnitude range. 

Figure [TT1 shows £ gm for "subhalo satellites", i.e., 
satellite galaxies hosted by intact subhalos (top panel) 
and for orphan satellites (bottom panel). For both 
types of satellite, the cross-correlation function is pos- 



Figure 10. Top panel: Distribution of host halo mass M200 
for central galaxies with —21 < M r < —20. Bottom panel: 
Mean host halo M200 (open diamonds) with 20-and 80- 
percentile values (dotted lines) are compared with the best-fit 
values recovered by modelling galaxy-mass cross-correlations 
(solid circles). 



itively biased with respect to £ mm on large scales, r ;> 
1 h~ 1 Mpc. This reflects the fact that relatively bright 
satellite galaxies almost all reside in halos more massive 
than M*, the characteristic mass at which halos are as 
strongly clustered as the underlying mass density field. 

On small scales, r <C 0.1 /i -1 Mpc, the shape of £ gm 
for satellite subhalos reflects the mass associated with 
the individual subhalos. A comparison between the £ gm 
curves in the top panel of Figure [11] and those in Fig- 
ure [9] shows that the average density profile of a satellite 
galaxy with a given luminosity is quite similar in both 
shape and amplitude to that of a central galaxy with the 
same luminosity. In detail, the subhalo satellite £ gm is 
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denser by ~ 10 — 20% with respect to the corresponding 
central galaxy £ gm , reflecting the hi gher concentration o f 
subhalos compared to main halos (IBullock et alj 12001). 

In the case of orphan satellites, negligible mass is 
associated with the satellite itself. On small scales the 
shape of the cross-correlation function reflects that of 
the host halos after convolution with the mean radial 
distribution of orphans within their hosts. Since or- 
phans are created by tidal stripping, they are found 
predominantly close to the centres of their hosts where 
tidal forces are strongest. Their r adial profile is t here- 
fore highly centrally concentrated (|Gao et al.ll2004l ) . Ta- 
ble [2] shows that the mean distance of orphan satellites 
from the centres of their hosts is significantly smaller 
than that of subhalo satellites. As a result, the orphan 
satellite £ gm is overdense on small scales with respect to 
fmm, since the one-halo term of the latte r reflects the 
halo density profile convolved with itself l|Ma and Fry! 

The cross-correlation with mass for all satellite 
galaxies is shown in Figure 1121 This is simply a linear 
combination of £ gm for the subhalo and orphan satel- 
lites, weighted by the relative fractions of each type. 
The shape of the resulting function appears quite sim- 
ple: unlike £ gm for the subhalo satellites, £ gm for the 
total satellite sample is positively biased with respect 
to £ mm on all scales, even for faint galaxy samples. In 
fact, its shape follows that of £ mm quite closely, with 
the addition of an upturn on small scales due to the 
mass associated with the individual subhalos. Encour- 
aged by this, we propose the following model for cross- 
correlations between satellite galaxies and the mass: 



£gm,sat (f) — 



Phalo(»";C, M) 



&(M host )£ mm (r) 



1 



3 exp 



(20) 



On scales r 3> 773, £ sa t is equal to the product of 
the mass autocorrelation function and the bias factor 
6(Mhost), where Mh os t is the average mass of halos which 
host satellite galaxies of a given luminosity. On scales 
r <C 2 /i _1 Mpc, a scale-dependent bias is apparent in 
the satellite £ gm , which we attempt to model with an 
exponential function. In eq. 1201 the parameters rp and 
3 control the characteristic scale and amplitude of this 
scale-dependent bias. On small scales, r < 0.1 /i _1 Mpc, 
the satellite model is dominated by the first term due 
to the subhalo density profile. 



Figure [12] shows the deviations of this model from 
the satellite £ gm measured in the simulation. Since a 
nonlinear model for £ mm is used as the basis of the satel- 
lite model on large scales, the systematic residuals in 
the fits to the halo and central galaxy cross-correlation 
functions are not present in the satellite model fits. The 
main deficiency in the satellite model is in reproducing 
the scale-dependent bias. The fits shown in Figure [T2] 
were obtained using single values of rp — 2 h~ 1 Mpc and 
3 — 0.5. They could be improved somewhat by varying 
these values for the different satellite galaxy luminosity 
bins, but systematic residuals still remain because the 
exponential function is only a crude match to the shape 
of the scale-dependent bias. 

The distribution of host halo masses for satellite 
galaxies in the luminosity range —21 < M r < —20 is 
shown in Figure 1131 The distribution spans five orders 
of magnitude in mass. Nevertheless, the best-fit model 
recovers the mean halo mass to within about 30%. The 
bottom panel shows mean and best-fit halo mass values 
as a function of M r . The values obtained from fitting 
the cross-correlations are accurate for host halo masses 
in the range —22 < M r < —19, but can differ from the 
true values by as much as 50% for host halos outside 
this intermediate range. 

The model parameters entering the first term in 
eq. [20] are the concentration, c, and the mass, M200 , 
of the halo density profile fitting formula. However, the 
best fitting parameter values (shown in Figure 1131 the 
concentration is taken as the typical value for halos of 
each mass) are not easily interpreted. 

4.3 Cross-correlations with mass for all 
galaxies 

In Figure [14] we present cross-correlations of all galax- 
ies with mass. This is simply the linear combination of 
the cross-correlations for central and satellite galaxies, 
weighted by the relative fractions of each type, i.e.: 



Cgm.all — (1 — /sat)£g 



ral ~\~ /sat^g 



(21) 



where the fraction of satellite galaxies is defined as 
/sat = A r sat/(A r ccntrai + AT sat ). The satellite fraction 
in each luminosity range strongly affects the shape of 
£ gm for galaxies of that luminosity. Table [2] shows that 
/sat ranges from ~ 40% in the lowest luminosity bins, 
M r > -19, to ~ 5% for M r < -22. 

Since eq.[2l]is a linear combination of two functions 
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Figure 11. Cross-correlation between the mass and satellite 
galaxies with (top panel) and without (bottom panel) associ- 
ated subhalos. In both cases, £ gm is positively biased on large 
scales, r > 1 h~ 1 yipc indicating that satellites preferentially 
reside in high mass halos. On small scales £ gm for subhalo 
satellites shows an upturn due to the mass bound to the in- 
dividual subhalos, whereas the orphan satellite £ gm reflects 
the density profile of the host halo convolved with the radial 
distribution of the orphan satellites. 



Figure 12. Top panel: Cross-correlations of the mass with 
all satellite galaxies. In this case the shape of § gm is similar 
to that of §mm with an upturn on small scales due to the 
subhalo mass associated with individual satellites that have 
not been fully tidally disrupted. Bottom panel: Deviations 
between the satellite £ gm and the model given by eq.[20] The 
deviations are largest at intermediate scales, r ~ 1 /i _1 Mpc, 
where the model fails to accurately describe the shape of the 
scale-dependent bias. 



which can differ in value by several orders of magnitude 
and since / sa t is of order unity, £ gm ,aii is dominated by 
the larger of £ gm , central and £ gm , S at at most radii. For in- 
stance, for —24 < M r < — 23, fgm,ceirtral/£gm,8at » 1 at 
all radii, thus eq. [21] gives £ gm ,aii ~ (1 - / sa t)£gm,ccntr a i 
and the cross-correlation function is very similar to 

£gm, central- At lower luminosities, £gm,sat/£ gm, central > 1 

at intermediate radii, 0.1 h~ 1 Mpc < r < 1 h~ Mpc, so 
£gm,aii ~ /satCgm.sat and in this case the cross-correlation 
function resembles that of the satellite galaxies. 

The combined £ gm model for central and satel- 



lite galaxies contains the following five parameters: 
Mcontrai, M host , M sat , c sa t , and / sat . However, in fit- 
ting this composite model to the measured cross- 
correlations, we obtained best-fit parameter values 
which were not in good agreement with the true val- 
ues of Mcentrai, Mhost and / sat . This suggests that the 
shape of £ gm is degenerate, with various combinations 
of the five free parameters giving almost equivalent fits. 

It is interesting to note that our cross-correlation 
model is based on components which deviate strongly 
from power laws, i.e., £n n , £mm and phaio- Nevertheless, 
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Figure 14. Top panel: Cross-correlations between all galax- 
ies and the mass at z = 0. High-luminosity galaxy samples 
are dominated by central galaxies, therefore £ gm resembles 
our simple model for £hm- Low-luminosity galaxy samples 
contain a significant fraction of satellite galaxies and £ gm fol- 
lows the shape of £ mm . Bottom panel: AT,(R) corresponding 
to f™. 



Figure 13. Top panel: Distribution of host halo M200 f° r 
satellite galaxies with —21 < M r < —20. Bottom panel: 
Mean host halo M200 (open diamonds) with 20-and 80- 
percentile values (dotted lines) are compared with best- 
fit values recovered by modelling the satellite-mass cross- 
correlations (solid circles). 



at intermediate luminosities £gm,aii is reasonably well 
described by a power-law over a wide range in radius. 
In order to illustrate this point, we plot in Figure [I5l the 
deviations from power-law fits to £ gm ,aii- In particular, 
we find that in the range —21 < M r < —20, £gm,aii is 
well fit by a power law with slope r -1,8 and the devia- 
tions from the fit are <C 10% at all radii! 

For galaxies more luminous than M r < —21, the 
deviations from best-fit power laws become J> 50% at 
some radii. Galaxies in this luminosity range are domi- 
nated by central galaxies, so our model for £ gm , central can 
be used to fit observations and provide detailed checks 
on various aspects of the cosmological and galaxy for- 



mation models. At lower luminosities, however, we con- 
clude that some independent information regarding the 
satellite fraction must be included in order to extract 
useful information about the average mass of the halos 
which host satellites. Alternatively, shear data can be 
stacked around galaxies that are observed to be brighter 
than all their neighbours, and thus are very likely cen- 
tral galaxies. Our models can be used to predict AE(J?) 
for direct comparison with such galaxy-galaxy lensing 
measurements. Examples are shown in the lower panel 
of Figure 1141 It is notable that these curves show a va- 
riety of scales and features which should be directly ob- 
servable. 



5 SUMMARY 

We have calculated cross-correlations between halo cen- 
tres and the mass, and between galaxies and the mass 
in the Millennium Simulation of a ACDM cosmology. 
The shape of the halo-mass cross-correlation function 
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Figure 15. Deviations from the best-fitting power laws of 
our measured cross-correlations between all galaxies and the 
mass. 

£hm is well fit by a simple two-part model. On small 
scales, £hm is specified by the halo density profile which 
is accurately described by the Einasto model advocated 
by N04. On large scales, is specified by the mass 
autocorrelation functio n predicted by linea r theor y and 
the halo bias model of ISheth and Tormenl (|l999T ) . The 
deviations between £hm and the best-fit model are dom- 
inated by the quasi-linear distortion in £ mm , but are 
otherwise < 5%. 

The cross-correlation functions of central galaxies, 
£gm, central, are reasonably well fit by our model for £ h m- 
The best-fit models recover the mean halo masses of 
central galaxies to within 30%. The cross-correlations 
of satellite galaxies, £ gm , S at appear qualitatively differ- 
ent from those of halos or central galaxies. Their shape 
is similar to that of £m m , with an upturn at small scales 
due to the mass associated with the individual subhalos 
in which most satellites reside. A model for £ gm , S at based 
on these features reproduces the cross-correlation func- 
tions to within 10% and recovers the mean host halo 
mass to within 50%. 

The cross-correlation function of all galaxies, £ gm ,aii 
is simply a linear combination of £ gm . central and £ gm ,sat 
weighted by the relative fractions of central and satel- 
lite galaxies. For very luminous galaxies, the satellite 
fraction / sa t <j 10%, and £ gm ,aii is dominated by the 
contribution of £ gm ,ccntrai • At intermediate luminosi- 
ties, —22 < M r < —20, £ gm ,aii is reasonably well fit by 
a power law. At lower luminosities the cross-correlation 
is dominated by the satellite galaxy contribution. 

The conversion from the three-dimensional cross- 
correlations to the directly observable mean tangential 
shear (which is proportional to AS(7?), the difference 
between the mean enclosed surface density and the lo- 
cal surface density at each projected radius R) accen- 
tuates features in the cross-correlations. Galaxy-galaxy 
lensing surveys typically contain enough information to 



separate the lenses into likely satellites and likely cen- 
tral systems. As a result, the features seen in our pre- 
dictions should provide information on cosmological pa- 
rameters, tidal stripping processes, and the exact way in 
which galaxies trace the dark matter. If the features are 
clearly detected where they are expected, this will pro- 
vide a major challenge to theories which try to replace 
dark matter by a modification of Einsteinian gravity. 
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